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ABSTRACT 

We use deep images acquired with the Advanced Camera for Surveys (ACS) on board of the 
Hubble Space Telescope (HST) in the filters F555W and F814W to characterize the properties of 
NGC 376, a young star cluster located in the wing of the Small Magellanic Cloud (SMC). Using 
isochrone fitting we derive for NGC 376 an age of 28 ± 7 Myr, in good agreement with previous 
studies. The high spatial resolution ACS data allow us to determine the center of gravity of 
the cluster and to construct extended surface brightness and radial density profiles. Neither of 
these profiles can be fitted with a theoretical model, suggesting that the cluster is not in virial 
equilibrium. Considering the young age of the cluster, we speculate that the distortion of the 
radial profiles may be the result of the rapid gas dispersal that follows the initial phase of star 
formation. The cluster shows clear evidence of dynamical mass segregation. From the properties 
of the radial profiles and the present day mass function (PDMF) we conclude that NGC 376 
appears to have already lost nearly 90% of its initial stellar mass, probably as a consequence of 
the sudden gas dispersal that follows the early phase of star formation (SF). 

Subject headings: (galaxies:) Magellanic Clouds - galaxies: star clusters: individual (NGC 376) 



1. Introduction 

Star clusters are expected to form over a 
wide spectrum of masses that is later modi- 
fied by the selective destr u ction of low-mass 
clust e rs ([Fall fc Reesl Il977t iGnedin fc Ostriker 
1997t IVesperini 19971 Il998l) . As summarized in 
Fall fc ZhaneJ (|200lh . star clusters are in general 
weakly bound and can be easily disrupted by a va- 



1 Based on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by AURA, Inc., under NASA 
contract NAS5-26555. These observations are associated 
with program 10248. 

2 Space Telescope Science Institute, 3700 San Martin 
Drive, Baltimore, MD 21218, USA 

3 European Space Agency, Research and Scientific Sup- 
port Department, Baltimore MD, USA 

4 Istituto Nazionale di Astrofisica, Osservatorio Astro- 
nomico di Bologna, Via Ranzani 1, 1-40127, Bologna, Italy 

5 Department of Astronomy, University of Wisconsin, 
475 North Charter St., Madison, WI 53706, USA 

6 Dipartimento di Astronomia, Universita degli Studi di 
Bologna, via Ranzani 1, 1-40127 Bologna, Italy 



riety of mechanisms that operate on different time 
scales. On short (t < 10 7 yr) and intermediate 
(10 7 < t < few x 10 8 yr) time scales, stellar evo- 
lution removes mass from star clusters by a com- 
bination of stellar winds, supernovae and other 
ejecta. On longer time scales (t > few x 10 8 yr) 
dynamical processes such as internal two-body re- 
laxation, gravitational shocks, and dynamical fric- 
tion become the main causes of mass loss. Mass 
loss through stellar evolution, in particular, seems 
to play a key role in shaping the mass function of 
the star cluster systems, from the initial power- 
law, that characterizes the mass function of the 
star forming regions, to the bell-sh aped mass func- 
tion of the globular clu sters (e.g. Vesperini ,1998; 



Vesperini fc Z cpf 2003) 



It is now believed that the fast disappearance 
of extremely young (< 10 Myr) clusters, ofte n 
known as "infant mortality" (|Lada fc Ladall2003h , 
is due to the r apid gas dispersa l caused by stel- 



lar winds 



Whitworth 1979). low- mass stars 



outflows ( Matzner fc McKeel2000f) , and early- type 
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star supernova explosions (jEgglestonl 120061 ) . The 
gas expulsion decreases the binding energy of the 
initial stellar system and, as a consequence, stars 
may suddenly have velocitie s that are highe r than 
the l ocal escape ve locity ( Tutukovl Il978 ; Hills 
1980tlGoodwirill987h . The probability that a star 
cluster will survive the gas dispersal phase de- 
pends largely on the efficiency of SF, a property 
that is still poorly quantifie d (lElmegreenl 120071 
Price fc Batell2009h . Recently. ISmith et all (|201ll ) 
have studied the effects of gas expulsion on sub- 
structured clusters formed under non-equilibrium 
initial conditions. They find that the initial spa- 
tial and kinematic distributions of the stars are far 
more important for cluster survival than the SF 
efficiency. They suggest that variations in clus- 
ter initial conditions, rather than SF efficiencies, 
determine whether clusters survive gas expulsion 
or not. In view of these findings, in depth ob- 
servational studies of the characteristics of young 
resolved star clusters are clearly needed. Measure- 
ments such as their PDMF, the presence of mass 
segregation, the stellar concentration and density 
distribution can provide much needed constraints 
for understanding the evolution and longevity of 
star cluster systems. 

Because of its close proximity (60.6 Kpc; 
Hilditch et"al~ll2005l h the SMC is uniquely suited 
for detailed investigations of the stellar content of 
regions of SF and young star clusters. With HST 
it is possible to spatially resolve the densest star 
clusters and perform a quantitative and accurate 
census of the stell ar content down to the sub-solar 
mass regime (e.g. Sabbi et all 12008) . As part of 
a project devoted to studying the properties of 
young and massive star clusters in the SMC (P.I. 
A. Nota, GO-10248), in this paper we present 
an in-depth study of the stellar content of the, 
so far, poorly investigated star cluster NGC 376 
(aj2ooo = 01 ft 03 m 53.5 s , Sj 20 oo = -72°49'27.0") 
based on HST imaging. NGC 376 is one of the 
brightest, richest, and youngest (from 16 Myr, 
Chiosi et al. 2006, to 25 Myr, Piatti et al. 2007) 
clusters in the eastern extension of the SMC to- 
ward the LMC, known also as the SMC wing. 
Previous studies suggest that the surface bright- 
ness profile o f NGC 376 depart s from a standard 
EFF model (jElson et al.l Il987l) a nd that it may 



The paper is organized as follows: in Section 2 
we present the observations and the analysis of 
the data, while Section 3 is dedicated to the de- 
scription of the color-magnitude diagram (CMD) 
and of the stellar populations found in the region. 
In Section 4 we discuss the structural parameters 
of NGC 376 as derived from the surface brightness 
and the stellar density profiles. In Section 5 we de- 
rive the age of the cluster from isochrone fitting. 
We present the luminosity function (LF) and the 
PDMF in Section 6. The results are discussed in 
Section 7. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

2.1. The data 

We used the HST/ ACS Wide Field Channel 
(WFC) to acquire deep images of NGC 376. The 
data were obtained in September 2004 as a part 
of a program devoted to study the characteris- 
tics of the youngest star clusters in the SMC 
(GO-10248; P.I. A. Nota). We acquired four 
450 sec long exposures (total integration time = 
1800 sec) in the F555W filter (~V), and four 490 
sec long exposures (total integrati on time=1960 
sec) in the F814W filter (^1 - see ISirianni et al 



be m erging with another cluster ([Carvalho et al 
2008h . 



2005, for a description of the ACS filters). Both 
F555W and F814W images were acquired follow- 
ing a four pointing dither pattern designed to bet- 
ter remove hot pixels and fill the gap between the 
two ACS/ WFC 2048 x 2048 CCDs. Dithering also 
allowed us to better sample the point-spread func- 
tion (PSF) and, by averaging the flat-field errors 
and smoothing over the spatial variation of the 
detector response, to improve the photometric ac- 
curacy. To recover the photometric information 
for the brightest, and otherwise saturated, stars, 
we acquired two short exposures (3.0 sec) in both 
filters. A summary of the observations is given in 
Tabled! 

All the images were acquired with a gain of 
2e~ ADLP 1 , and the entire data set was processed 
through the standard Space Telescope Science In- 
stitute ACS calibration pipeline CALACS to sub- 
tract super-bias, and super-dark, and to apply 
the flat-field correction. Our ACS/WFC images 
cover an area of 200" x 200" and have a pixel 
scale of 0.05" pixel -1 . Assuming a distance mod- 
ulus (m — M)o = 18.912 ±0.035 (corresponding to 
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Table 1 

Journal of WFC/ACS observations 



Image Name 


Date & Time of Observation 


R.A. 


Dec. 


Filter 


Exposure Time 


TQ9FD8T, ZD 


12/09/04 


13:52:45 


ni hnqmcqso 




-T OOO VV 


o .u 


J92F08M0Q 


12/09/04 


13:55:06 


01 ft 03 m 53 i ?8 


-72°49'29"6 


F555W 


450.0 


J92F08M2Q 


12/09/04 


14:05:12 


01 h Q3 m 53!8 


-72°49'29'.'6 


F555W 


450.0 


J92F08M4Q 


12/09/04 


14:15:18 


01 h 03 m 53!8 


-72°49'29'/6 


F555W 


450.0 


J92F08M6Q 


12/09/04 


14:25:24 


01 h 03 m 53 i :8 


-72°49'29'/6 


F555W 


450.0 


J92F08M9Q 


12/09/04 


14:41:01 


01 h 03 m 53!8 


-72°49'29'/6 


F555W 


3.0 


J92F08MBQ 


12/09/04 


14:49:52 


01 h Q3 m 53!8 


-72°49'29'/6 


F814W 


3.0 


J92F08MCQ 


12/09/04 


14:52:13 


01 h 03 m 53!8 


-72°49'29"6 


F814W 


490.0 


J92F08MEQ 


12/09/04 


15:29:33 


01 h 03 m 53!8 


-72°49'29'/6 


F814W 


490.0 


J92F08MGQ 


12/09/04 


15:40:19 


01 h 03 m 53!8 


-72°49'29'/6 


F814W 


490.0 


J92F08MIQ 


12/09/04 


15:51:05 


01 h 03 m 53!8 


-72°49'29"6 


F814W 


490.0 


J92F08MLQ 


12/09/04 


16:07:22 


01 h 03 m 53!8 


™72°49'29"6 


F814W 


3.0 



- 60.6 kpc. Iffilditch et al.l l2005h the ACS/WFC 
field of view corresponds to a projected area of 

— 50 x 50 pc 2 . A color composite ACS/WFC im- 
age of NGC 376 is shown in Figure [1] 

2.2. Photometric Reduction 

The photometric analysis was carried out di- 
rectly on the pipeline-corrected images with extcn- 
sion _flt, using the prog ram img2xym_WFC. 09X10 
( Anderson k, Kinj2006 ). that was specifically de- 
signed to perform photometric analysis of the 
under-sampled ACS/WFC data. A library of em- 
pirical PSFs, that takes into account the PSF spa- 
tial variation due to the telescope o ptics and th e 
variable charge diffusion in the CCD (|Kristll2003h . 
is provided with the code. Temporal variations 
in the PSF, caused by changes in the telescope's 
focus induced by spacecraft breathing, are also 
taken into account by fitting in each image the 
library PSF to the brighter (S/N > 100) stars. 
We used the routine described in lAnderson et al 



(2008) ;o further refine the photometry. The pro- 
gram divides all images in to regions of 25 x 25 pix- 
els and, region after region, simultaneously finds 
and fits the stars in all the exposures. For each 
star found, the code provides the average X and Y 
coordinates and F555W and F814W magnitudes, 
the r.m.s. of coordinates and magnitudes, as well 
as an estimate of the psf-fitting quality. 




Fig. 1. — ACS/WFC color-composite image of 
NGC 376. The image acquired through the filter 
F555W is in blue, while the F814W image is in 
red. North and East directions are indicated. The 
solid red circle defines the area within the core ra- 
dius, the yellow solid circle marks the area within 
the tidal radius, and the green solid circle shows 
the extent of the cluster's tail. 
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To remove as many spurious detections as pos- 
sible from the final catalog, we retained only those 
stars that were identified in both the filters in 
at least three exposures with a positional error 
smaller than < O.f pixel. Photometric r.m.s as 
a function of magnitude are shown in Figure [2j 

To calibrate the final catalog in the Vegamag 
photometric system we combine d the _flt images 



using the multidrizzlc package (jKoekemoer et al 



2002) and performed aperture photometry for sev- 
eral isolated stars for both the F555W and F81 4W 
filters using DAOPHOT in IRAF0. These stars 
were calibrated into the Veg amag photomet r ic sys- 



(2005) 



tern following the recipe by ISirianni et al. 
and then used to calibrate the photometric cata- 
log. 

2.3. Artificial Star Tests: Photometric Er- 
rors and Completeness 

Artificial star tests are a standard procedure 
used to quantify the completeness of a photo- 
metric catalog, as well as to test the impact of 
crowding on the photometric accuracy. The tests 
are performed by inserting stars with known flux 
and position in the data set, and then repeat- 
ing the photometric analysis using the same pro- 
cedure applied to the real data. The difference 
between the input and output magnitude of the 
recovered artificial stars (Figure [3]) provides an 
estimate of the photometric accuracy, while the 
photometric completeness, as a function of mag- 
nitude, is derived from the fraction of recovered 
artificial stars with respect to the simulated ones. 
Figure @] shows the variation of the completeness 
factor [Cf ) as a function of the distance from the 
center of NGC 376. 

The same routines used in Section 12.21 to re- 
fine the photometry, can also perform artificial 
star tests. The program divides each image in 
to regions of 25 x 25 pixels, and adds an artifi- 
cial star in a region at a time. As during the 
photometric analysis, the program finds and fits 
all the sour ces in the region simu ltaneously in all 
the frames ( Anderson et al.ll2008l ). This approach 
avoids the problem that artificial stars interfere 
with each other, artificially increasing the crowd- 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by AURA, Inc., under co- 
operative agreement with the National Science Foundation. 
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Fig. 2. — Photometric r.m.s as a function of mag- 
nitude for stars that have been detected in at least 
three exposures after the selection in positional er- 
rors, in the F814W (upper panel) and the F555W 
(lower panel) filter respectively. 
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Fig. 3. — Difference between the input and output 
magnitude of the artificial stars as a function of 
the input magnitude for both the F8f4W (upper 
panel) and the F555W (lower panel) filters. 
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ing of the image, and thus altering the inferred 
completeness. In total we simulated more than 
1,000,000 stars in each of the F555W and F814W 
exposures. 

We considered an artificial star as recovered if: 



• The input and output fluxes agree to within 
0.75 mag; 

• The input and output positions agree to 
within 1 pixel; 

• The star is found in at least three F555W 
and three F814W exposures, with a posi- 
tional error < 0.1 pixel. 

We find that the completeness of our sample 
varies as a function of the distance from the clus- 
ter (Figure [4])- For example within the innermost 
5" from NGC 376 the sample is 50% complete to 
"T-F555W = 25.25, between 5 and 10" the same 
completeness is at mF555w — 25.5, and it goes 
down to mp555w = 26.25 between 10 and 15". At 
WF555W = 24.75 and mF8i4W = 24.25 the sample 
is > 65% complete over the entire field. 

3. COLOR-MAGNITUDE DIAGRAMS 

The CMD mp8i4W versus mpsssw — "t,F814W 
of all the stars that passed our selection criteria 
i s shown in Figure [5] At a distance of 60.6 kpc 



(Hild itch et al.l I2005T) and assuming a metallicit y 
Z=0.001 dCarrera et al.l 120081 IParisi et al.ll2010h . 
Padua isochrones (|Bertelli et al.l 120081 120091) pre- 
dict that the 13 Gyr main sequence turn-off 
(MSTO) is at m F8 i 4W = 21.88. Our photom- 
etry therefore reaches > 4 magnitudes below 
the expected oldest MSTO. Padua isochrones 



(|Bertelli et al.l 120081 120091 ) for Y=0.24 and dif- 
ferent ages and metallicities are superimposed 
on the CMD of Figure [5] for reference. To plot 
the isochrones we assumed a distance modulus of 
(to — M)o = 18.9, the Galactic extinction law, a 
reddening of E(B-V) = 0.08 and A F555W =0.19. 

An inspection by eye of Figure [5] shows that, 
as found in several other SMC stellar fields (i.e . 
Sirianni et al.l 120021: iMcCumber fc Garnettl [2 005: 



Nota et all2006llSabbi et alJ2007tlChiosi fc Vallenarl 
' Carlson et al.ll2007l : lNoel et al.l2007USabbi et al.l 



200 




Fig. 4. — The Completeness Factor Cf as a func- 
tion of the input magnitude as derived from the ar- 
tificial star tests for both the F814W (upper panel) 
and the F555W (lower panel) filters as a function 
of the distance from the NGC 376 center. The 
solid lines represent the CfS within 5" from the 
center, the Cps in the annulus between 5" and 
10" are represented by the dotted lines, the dashed 
lines are the CfS between 10" and 15", the long- 
dashed lines, the dashed-dotted lines and and the 
long-dashed-dotted lines are the CfS between 15" 
and 25", 25" and 35" and 45" and 55" respectively. 



2009), stellar populations of different ages coexist 



in the area. 
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Fig. 5.— ACS/WFC CMD 

"T-F814W VS. TOF555W — 

WF814W of NGC 376 with superimposed Padua 
isochrones of different ages and metallicities. The 
blue continuous isochrone corresponds to a 28 Myr 
old stellar population with metallicity Z=0.004, 
the dotted blue line corresponds to a Z=0.004 200 
Myr old stellar population, and the dashed blue 
isochrone is for a Z=0.004 500 Myr old population. 
The green continuous and the long-dashed lines 
correspond to Z=0.002 isochrones for a 3 and 5 
Gyr old stellar populations respectively. Finally 
the red line is a Z=0.001, 8 Gyr old isochrone. 



Old stars: A stellar population older than 
~ 2 Gyr can be easily distinguished in the CMD. 
Evolutionary phases associated with this popula- 
tion are: 

• The well defined lower main sequence (MS), 
which extends from m.F8i4W — 21.7 down to 
m-F8i4W < 26.4; 

• The broad (20.3 < mpsuw 21.6) subgiant 
branch (SGB), visible in the color range 
0.42 < m F5 5 5W — m,F8i4W 0.90. Compar- 
ison with the isochrones suggests that the 
majority of the stars in this evolutionary 
phase are ~ 4 — 5 Gyr old; 

• The bright red giant branch (RGB), with 
the brightest star at mF8i4W — 14-9 and 
WF555W — WF814W — 1-94. The relatively 
small number of RGB stars does not allow 
us to establish if the brightest stars of the 
RGB correspond or not to the RGB-tip; 

• The red clump (RC), at TOF814W — 18.5. 

Asymptotic giant branch (AGB) stars are likely 
present in the CMD, but it is difficult to separate 
them from RGB stars. 

Young and intermediate-age stars: The 
bright (15.0 < mp8i4w S 21.0) and blue (mF555w — 
WF814W < 0.2) well-populated upper MS indicates 
that in this region the SMC was still forming stars 
less than ~ 2 Gyr ago. Young low mass stars 
are likely present in the lower MS, however they 
cannot be distinguished from the older (> 2 Gyr) 
MS. This finding is also supported by the pres- 
ence of stars, above the RC, in the magnitude 
range 16.5 < ?TiF8i4W 18-0 and color range 
0.56 < m.F555w — "IF814W J$ 0-89 which likely cor- 
responds to the lower end of the blue edge of the 
blue loop and thus are younger than ~ 500 Myr. 

Between 15.0 < m F8 i4w < 17.3 and 0.0 < 
"T-F555W — "^F8i4W ;$ 0.16 there is a secondary 
sequence of objects that are likely Be stars. Their 
presence in the CMD suggests that this region 
was s till forming stars ~ 30 Myr ago (jKeller et al. 
2000). This is also confirmed by the few red super 



giant stars (RSGs) brighter than ttifsmw > 13.1 
and redder than mF555w — "iF8i4w > 1-24- 

In two other SMC star clusters (namely NGC 34 6 
and NGC 602, iNota et al-EooSlSabbi et al.ll2007 : 
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Carlson et all l2007t ICignoni et all 120091 |201lh a 

similar observational set-up allowed us to find 
rich populations of pre-main sequence stars on 
the right side of the lower MS. The absence of 
pre-main sequence stars in Figure [5] suggests that 
in this region of the SMC very few stars, if any, 
were formed in the last ~ 15 Myr. 

Finally we do not find any clear evidence of 
horizontal branch (HB) stars in the CMD of Fig- 
ure which are considered unequivocal indicators 
of a metal poor stellar population, older than ~ 
10 Gyr, however the presence of an old-metal poor 
population in the SMC is confirmed by the find- 
ing of several RR-Ly rae over the entire SMC (e.g. 
ISosvnski et al. 2002 ). The paucity of HB stars 
has been confirmed by other photometric stud- 
Noil et alJl2007t ISabbi et al 



i.e. 



iesof the SMC 
2009), and this apparent dichotomy between the 



presence of RR-Lyrae and the absence of extended 
HB has been interpreted by several autho rs (i.e. 



Chiosi fc Vallenarill2Q07t iNoel et all lioioh 



as an 

indication that the SMC formed very few stars in 
the first 2-3 Gyr . 

4. NGC 376 STRUCTURAL PARAME- 
TERS 

Once formed, the evolution of a star cluster is 
affected by a continuous mass loss caused by gas 
expulsion, low-mass star evaporation and stellar 
evolution. Surface brightness and stellar density 
profiles are commonly used to probe the dynami- 
cal status of a star cluster. The first step in build- 
ing a radial profile is to determine the center of 
the stellar population. Since stellar luminosities 
are not always proportional to stellar masses, the 
center of luminosity (C/ um ) can differ significantly 
from the center of gravity {C grav ) of a cluster. 

In the F555W image the Ci um is at a = 
01 h 03 m 53.9 s ± 0.01 s ; 6 = -72°49'34.0" ± 0.3". 
We exploited the high spatial resolution of our 
photometric catalog to determine also the C grav 

we itera- 



of NGC 376. 



To find the C 



tively averaged the stel lar X and Y coordinates 



(|Montegriffo et al.l Il995h . In order to take into 



account effects due to crowding and incomplete- 
ness, we computed the C grav using three dif- 
ferent magnitude thresholds (mF555w < 22, 23, 
and 24). Finally to test the impact of the SMC 
field on our measurement, we derived the C gra v 



both using all the stars, and using only stars 
bluer than ttifsssw — "T-F814W < 0.5. Our fi- 
nal best estimate of the NGC 376 C grav is a — 
0l' l 03 m 52.98 s ± 0.01; 5 = -72°49'32.5" ± 0.15". 



4.1. The surface brightness profile 

NGC 376 has been recently observed in the V 
band wit h the ESQ Danish 1 .54 m telescope, in 
La Silla (jCarvalho et alj|2008j) . The derived sur- 
face brightness profile shows a central peak, and 
an external "bump" that can not be reproduce d 



with a standard EFF model (Els on et alJ ll987) 



Carvalho et al. (2008) interpreted these anomalies 
as the result of a recent, or even still ongoing, 
merger with another star cluster. 

To measure the surface brightness profile of 
NGC 376 (Figure [6]) we divided our catalog into 
equally spaced annuli centered on the C; um . Be- 
cause of the cluster location in the upper right 
corner of the image (Figure [!} only the annuli 
in the range 0" — 50" are fully imaged, while at 
larger radii only a portion of each annulus falls in 
our field of view. For this reason we divided each 
annulus into 12 sectors, and for each annulus we 
considered only those sectors that fully lie in our 
field of view. 

Figure |4] shows that for stars brighter than 
WF555W ^ 24 our catalog is more than 60% 
complete at each distance from the cluster cen- 
ter. We, therefore, used only stars brighter than 
this threshold to measure the surface brightness 
profile. Figure 0] shows also that, because of 
the increasing crowding, the completeness correc- 
tion changes as a function of the distance from 
NGC 376. To apply the right correction to each 
point of the profile, in each annulus we divided the 
stars of each sector into bins of magnitude, and to 
each magnitude bin we applied the completeness 
factor as measured at that distance from the 
cluster center. 

We used the recipe from ISirianni et al. (2005) 
to convert the F555W magnitudes into fluxes. The 
brightness of each sector was obtained by summing 
the flux of all the stars brighter than ttifsssw < 24. 
The flux of each annulus was obtained from the 
average of its sectors. 

The well defined RGB and SGB in Figure [5] in- 
dicate that a considerable number of stars in our 
catalog belong to the field of the SMC. To statisti- 
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Fig. 6. — NGC 376 surface brightness profile with 
the EFF models superimposed. The black line rep- 
resents the best fit to our data, while the gray 
region represents the uncertainties to the fit. In 
the Top Panel the average surface brightness of 
the field of the SMC has been measured from 
SF10. Red dots are for all the stars brighter than 
WF555W < 24, while blue triangles are derived con- 
sidering only stars between 14 < tofsssw < 24. 
In the Bottom Panel the surface brightness of the 
SMC field has been estimated from the stars at a 
distance larger than 120" from the cluster. 



cally remove the contribution of the SMC from the 
surface brightness profile, we measured the aver- 
age surface brightness in SF10, a field in the wing 
of the SMC (a = 0l' l 08 m 39 s ; S = -72°58'45") 
that was observed in January 2006 with a simi- 
lar observational strategy with the ACS/WFC in 
the F555W and F814W filters (GO-10396, P.I. J.S. 
Gallagher) , as part of a project devoted to study- 
ing the star formation history (SFH) of the field 
of the SMC (|Sabbi et al.ll2009h as well as to char- 
acterizing th e properties of inte rmediate and old 
star clusters (jGlatt et alJl2008al |b[). 

The surface brightness profile of NGC 376, after 
the subtraction of the stars SMC field, is shown in 
the top panel of Figure |5] (red dots). To fit this 
profile we used an EFF model expressed in mag- 
nitude per surface area, rather than luminosity, 
using the formula 

lfi(r) = lfj, (r) + 1.257iog(l + r 2 /a 2 ) 

where Ifio is the central surface brightness (corre- 
sponding to —2.5Loqiio(r) in t he original formula 
presented in lElson et al.lll987l ). 7 is the slope of 
the power-law, and a is a dimensionless param- 
eter that is related to the core radius r c of the 
King profile by r c = a(2 2 ^ - l) 1/2 (= 8.10 ± 1.58 
arcsec). The surface brightness profile is quite ir- 
regular, and shows an evident central peak that 
cannot be reproduced by our best fit. In agree - 
ment with the findings of ICarvalho et al.1 (|2008l ). 
beyond ~ 60" there is a bump in the counts that 
exceeds the EFF profile. 

To verify that the irregularities found in the sur- 
face brightness profile are not caused by few bright 
foreground stars nor by local variations in the field 
of the SMC we repeated the analysis by first se- 
lecting only the stars fainter than ttifsssw > 14.0 
(blue triangles in the top panel of Figure [6]), and 
then by measuring the average surface brightness 
at a distance larger than 120" from NGC 376 Ci um 
(Figure [6]- bottom panel). 



4.2. The stellar density profile 

Surface density profiles are dominated by the 
light from the brightest stars, but low mass stars 
play a dominat role in the dynamics of a star clus- 
ter. To better constrain NGC 376 dynamics, and 
taking advantage of the high spatial resolution of 
the ACS data, we also analyzed the stellar density 
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profile of the cluster. Following our adopted pro- 
cedure for the surface brightness profile, the stellar 
density profile was obtained by dividing our sam- 
ple in equally spaced annuli centered, in this case, 
on the Cg r av Each annulus was divided into 12 
sectors to take into account the fact that, at a dis- 
tance larger than r > 50", part of the annuli fall 
outside the image. In our analysis we considered 
only those sectors that fully lie in our field of view. 
In each sector we divided the stars brighter than 
WF555W ^ 24 into bins of magnitude and then we 
applied the completeness factor Cf, as measured 
at that distance from the cluster center, to each 
magnitude bin. We obtained the stellar density of 
each sector by dividing the number of stars by the 
area of the sector. The average stellar density of 
the sectors in an annulus was used as the stellar 
density of that annulus. To estimate the density in 
the field of SMC we used both the average stellar 
density of SF10 (Figure[7]- Top Panel) and the av- 
erage stellar density at a distance larger than 120" 
from the cluster C grav (Figure 0- Bottom Panel). 
The radial density profile of NGC 376, after the 
subtraction of the SMC background, is shown in 
Figure E 

King models ( Kind 19661) can be used to infer 
the dynamical status of a star cluster. King pro- 
files are usually described by a core radius r c , that 
corresponds to the distance from the center where 
the projected stellar density falls to 0.5013, and 
a tidal radius r t , that indicates where the poten- 
tial well of the cluster equals zero. The ratio of 
the two radii (r t to r c ) defines the concentration 
of the star cluster (c = log 10 (r t /r c )) and it can 
be used to establis h if a stellar system is in virial 
equilibrium or not dMevlan fc Heggielll997t ): 

• Models with c between 0.75 and 1.75 fit re- 
laxed star clusters very well. 

• Models with c > 2.2 fit core collapsed glob- 
ular clusters. 

• Stellar systems described by a King model 
with c < 0.7 are not in virial equilibrium. 

The King model that best fits our data (r c = 
7.95 ± 0.35; r t = 20.0 ± 0.67) provides a quite 
low concentration (c ~ 0.4), suggesting that the 
cluster is not in virial equilibrium. To facilitate 
a more direct comparison between the EFF and 




0.5 1 1.5 2 
Log(r}[arcsec] 




0.5 1 1.5 2 
Log(r)[arcsec] 



Fig. 7. — NGC 376 stellar density profile with su- 
perimposed the King profile (black line) in better 
agreement with the data. The grey region repre- 
sents the uncertainty of the fit. The local aver- 
age SMC stellar density has been estimated both 
from SF10 (Top Panel) and at a distance larger 
than 120" from the cluster (Bottom Panel). The 
dashed blue line in the bottom panel is EFF model 
that best fit the surface density profile shown in 
Fig. |6l scaled to the stellar density. The inconsis- 
tency of the density profile with both the King's 
and EFF models is apparent. 
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King models, in Figure [7] (Bottom Panel) we plot- 
ted also the best fit EFF model, derived in Sec- 
tion 14.11 This comparison shows that a single 
King (or EFF) model cannot fit the entire stel- 
lar density profile. Similar to what we found for 
the surface brightness, our best fit underestimates 
the stellar density in the innermost ^3 — 4", and 
cannot reproduce the tail of counts that contin- 
ues well beyond (up to ~ 60") the tidal radius 
(r t = 20 ± 0.67") estimated from the best fit of 
the stellar density profile. 

Dynamical simulations ( Kupper et al.ll201o( ) us- 
ing the velocity dispersion profile show that the 
number of potential escapers from the star clus- 
ter increases and, sometimes even dominates the 
surface density profile, for radii larger than ~ half 
the Jacobi radius (often approximated to the r t of 
the King model). In this context the comparison 
between the surface brightness and the stellar den- 
sity profiles seems to suggest that a considerable 
fraction of the stars that were formed in NGC 376 
are not gravitationally bound to the cluster any 
more. 

5. THE AGE 

Figure [5] -panel (A) shows the CMD of all the 
stars brighter than mF8i4W < 23 found within the 
NGC 376 tidal radius. The most evident feature 
in the CMD is the tight blue (mF555w — WF814W — 
0.0) MS, likely populated by NGC 376 stars. How- 
ever the presence of SGB and RGB stars redder 
than mFsssw — wfsmw > 0.4 indicates that even 
inside the cluster the contribution of stars belong- 
ing to the field of the SMC is not negligible. To 
characterize the field of the SMC around NGC 376 
we selected stars brighter than tofsi4W < 23 that 
are at a distance r > 120" from the C grav of the 
cluster. This selection allowed us to avoid the pop- 
ulation of "potential escapers" that likely domi- 
nates between - 20 - 60" (Figured]). The CMD 
of the field, normalized to the area covered by 
NGC 376, is shown in Figure [5]- panel (B), while 
the CMD of NGC 376 after the subtraction of the 
SMC field is shown in Figure H- panel (C). 

W e used Padua isochrones ( Bertelli et al.ll2008l 
2009) computed for Y=0.24 and various Z to in- 
fer the age of NGC 376 from the mF8i4W vs- 
WF555W — "1F814W CMD. To reproduce the tight 
MS of NGC 376 we started from the literature 



(A) 




NGC 376 - SMC field 



(C) 



0.5 1 1.5 



0.5 1 1.5 
^F555W ~~ m F814W 



0.5 1 1.5 



Fig. 8. — CMDs mF814W VS. TOF555W ~ TIF814W of 

the stars within the tidal radius of NGC 376 (panel 
(A)), the SMC field, as observed at r > 120" from 
the center on NGC 376 and normalized to the area 
of the cluster (panel (BJ). Panel (C) shows the 
CMD of NGC 376 after the decontamination from 
the SMC field. 
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values for distance, reddening and metallicity, and 
then we varied one parameter at the time. The 
effects of these changes are shown Figure El 

We chose for the di stance modulu s of the SMC 
the value derived by iHilditch et all (|2005t ) from 
the analysis of more than 50 eclipsing OB star 
binary systems (to - M) = 18.912 ± 0.035. 
The binaries are scattered over the whole galaxy, 
and therefore are likely representative of the 
galaxy mean distance. It has to be noted 



that several authors ( e.g. |M athews on et al 



Hatzidimitriou et al 



1988; 



2001 



. _ - . 1993: ICrowl. et al . _ 

Lah. et alJ2005HGlatt et al.l2008al) found that the 
depth of the SMC can be up to 20 kpc. Lithe sam- 
ple of clusters analyzed bv ldatt et al. ( 2008a ). for 



example, the closest object has a distance modu- 
lus (to — M)o = 18.5, while the distance modulus 
of the farthest system is (m — M) = 19.12. Panel 
(A) of Figure[5]shows that among this range of val- 
ues our data are well fitted by the average distance 
modulus (to — M)o — 18.9. Similarly panel (B) of 
Figure [9] shows that the average SMC r eddening 
value E(B-V)=0.08 (jZaritkv et i alJl2002h matches 
our data well. We also considered isochrones 
for three different metallicity values (Z=0.001, 
Z=0.004 and Z=0.008 - Panel (C) of Figure© and 
found that a metallicity of 1/5 solar is the most 
adequate to reproduce both the colors and mag- 
nitudes of NGC 376 MS stars from tofsi4W ~ 16 
down to TOF814W ~ 23. In summary, assuming 
a distance modulus (to — M)q = 18.9, a redden- 
ing E(B - V) = 0.08 and a metallicity Z=0.004, 
Padua isochrones indicate that NGC 376 has an 
age of 28 ± 7 Myr (Figure 0- panel (D)), i n good 
agreement with the results by lPiatti et al.l ( 2007 ) 
and with the finding that NGC 376 is hosting 
several Be stars. 



6. PRESENT DAY LUMINOSITY AND 
MASS FUNCTIONS 

We now examine the stellar PDMF of NGC 376 
and derive the total mass of the cluster. We used a 
procedure analogous to that described in Section|4] 
to measure the LF of NGC 376 down to TOF555W = 
24 and to apply a correction for completeness as 
a function of the distance from the cluster. To 
remove the contribution of the SMC stars from 
the LF we used the LF of the stars found at a 
distance larger than r > 120". The catalog was 




o 0.5 1 1.5 



0.5 1 1.5 



"FSBBW '"fbuw 



Fig. 9.— CMD TOF814W VS. m F5 55W ~ ?™F814W 

of NGC 376, after the subtraction of the SMC 
field, with a 28 Myr Padua isochrone computed for 
Z=0.004, assuming an E(B-V)= 0.06 and differ- 
ent distance moduli (panel (A)); a 28 Myr Padua 
isochrone computed for Z=0.004, assuming a dis- 
tance modulus (m-M)o = 18.9 but different values 
of reddening (panel (B); 28 Myr Padua isochrones 
computed for different metallicities, assuming a 
distance modulus (m-M) = 18.9 and an E(B- 
V)= 0.06 (panel (C)). Panel (D) shows Z=0.004 
Padua isochrones in the interval of time between 
32.4 and 35.5 Myr assuming a distance modulus 
(m-M)o = 18.9 and an E(B-V)= 0.06. 
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corrected for completeness and normalized for the 
area. The LFs in the F555W band as a function 
of the distance from NGC 376 C gra v are shown in 
Figure [TU] (left column) before (black histogram) 
and after (red histogram) the subtraction of the 
SMC background. 

To convert the observed LFs in to PDMFs we 
used the mass-luminosity (ML) relation derived 
from the 2 8 Myr old, Z=0.00 4 Y= 0.24 Padua 
isochrone (|Bertelli et al.l 120081 l2009h . Among 
the available parametrizations, we chose the one 
proposed by IScalol ( 19861 ). in which the MF 
is characterized by the logarithmic derivative 
T = dlog£(log(ra))/d, where £(log(ra)) is the MF 
and r is its slope. In this parametrization, the 
slope of the sol ar neighborhood initial mass func- 
tion derived by ISafpeterl (Il955h is r = —1.35. 

We calculated the derivative of the ML rela- 
tion in the mass range between 1-10 M Q using 
a spline interpolation. The PDMFs derived in 
this way for different radii are shown in Figure [10] 
(right column). The derived slopes are affected 
by several uncertainties, such as the assumptions 
on distance, amount of extinction, and age as 
well as residual contamination from the field. A 
major source of uncertainty for NGC 376 likely 
comes from the unresolved binary systems. In 
a dynamically mass segregated cluster, the most 
massive binaries should be confined in the cen- 
ter with respect to light er systems, further steep- 



ening the distribution. ISagar fc Richtlerl ((.1991) 



estimated that if each star in the mass range 
2-14 Mq has one companion, the average IMF 
slope derived for five young clusters in the LMC 
would si gnificantly steepen. More detailed mod - 
eling bv lWeidner. Kroupa fc Maschbergen (|2009h . 
however, suggests that binaries have little effect 
on the measured slope of the MF. 

A weighted least-mean square fit of the data 
indicates that the PDMF remains flat ((r) = 
—0.31 ± 0.27) over the entire cluster extent, and 
it becomes close to the value found by Salpeter 
only in the c luster's tail. As for other young 
star c lusters (Ide Griis et al.l I2002E ISirianni et al 



20021: IStolte et al.l l2002t iGouliermis et al.l 120041 
Sabbi et al.l 120081 ) the steepening of the PDMF 
with distance from the center is caused by a lack 
of massive stars rather than an excess of low- 
mass stars in the outskirts. The paucity of mas- 
sive stars in the halo of a star cluster is normally 
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Fig. 10. — LFs (left column) and MFs (right col- 
umn) for different annuli (0"-10"; 10" -20"; 20"- 
30"; 30" - 40"; 40" - 50") around NGC 376 C grav . 
In the left column black histograms are the LFs, 
corrected for completeness, before the SMC back- 
ground subtraction, and red histograms are the 
LFs after the subtraction of the SMC field. The 
slope T of the MF, as derived from the weighted 
least squares fit to the data, is shown in each of 
the right panels. 
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interpreted as a signature of mass segregation. 
Some of the known mass-segregated clusters are 
so yo ung that mass-seg r egation has to be primor- 



dial ([Hillenbrand 1997]; iHillenbrand k Hartmann 
19981; ISabbi et al.ll2008h . This does not seam to be 
the case for NGC 376: from the PDMF and the 
King profile we find that the stellar mass within 
the tidal radius r t is MroT rt — 7 x 10 4 M Q , and, 
considering also the stars in the cluster's tail (r < 
60"), we find a total mass M TOT ~ 6 x 10 5 M . It 
thus appears that only ~ 10% of the present day 
cluster mass is within the tidal radius. We can 
substitute the cluster virial radius with the half 
mass radius (r V i r r/ lm ), in which case (consider- 
ing only the stars within the tidal radius) the dy- 
namical time scale (Td yn ) necessary for NGC 376 
to reach dynamical equilibrium is Td yn ~ 2.6 x 10 5 
yr , and, if we consider the mass within 60" from 
the Cg r av, it rises to Ta yn — 7.8 x 10 5 yr , pointing 
toward a dynamical origin for the mass segregation 
observed in NGC 376. The PDMF also suggests 
that even within the tidal radius stars are weekly 
bound since at the tidal radius (which at the dis- 
tance of the SMC correspond to ~ 5.9 pc) the 
escape velocity is ~ 3km/s. These findings sug- 
gest that the NGC 376 is likely evaporating into 
the field of the SMC. The properties of NGC 376 
derived in this paper are summarized in Table [2j 

7. DISCUSSION AND CONCLUSIONS 

We analyzed deep ACS/WFC images of the 
young star cluster NGC 376 in the filters F555W 
and F814W, as part of a project devoted to study 
the formation and evolution of young star clus- 
ters in the SMC (P.I. A. Nota; GO-10248) Us- 
ing Padua isochrones (jBertelli et all 120081 l2009h 
we found that NGC 376 is 28 ± 7 Myr old, in 
agreem ent with the previous study bv lPiatti et al 
(|2007t ). 



Previous analysis of the cluster surface bright- 
ness profile (jCarvalho et al.l 120081 ) showed that it 
is quite irregular, with a spike in the center and 
an extended bump in the outer regions. We took 
advantage of the high spatial resolution of our 
dataset to repeat the analysis of the surface bright- 
ness profile and to measure, for the first time, 
the cluster stellar density distribution down to 
WF555W = 24, which, at the distance of the SMC, 
for a 28 Myr old stellar population of Z=0.004 



metallicity, corr esponds to c± 0.96 M Q . We used 
an EFF model (jElson et al.l Il987h to fit the s ur- 
face brightness and a King model (|Kinglll966l ) to 
reproduce the stellar density profile. In both the 
cases, the best fits underestimate the central peak. 
Furthermore the stellar density profile shows a tail 
(r ta ii — 60" ~ 17 pc) that extends well beyond the 
tidal radius (r t = 20" ~ 5.9 pc) derived from the 
best fit King model. Our analysis indicates that 
NGC 376 has a low concentration (c = 0.4), typ- 
i cal of clusters that are not in virial equilibrium 
(|Mevlan k Heggielll997h . 

The anomalous surface bright ness profile of 
NGC 376 was interpreted by ICarvalho et al 



(|2008l ) as the result of a recent, or even a still 
ongoing, merger. However, no obvious interacting 
candidate, such as another star cluster was found 
in the vicinity of NGC 376, and the tight MS of 
NGC 376 suggests that all its stars were formed in 
a relatively short interval of time. Similarly we did 
not find any know giant molecular cloud (GMC) 
within ~ 30', which at the distance of the SMC 
corresponds to a pr ojected distance of ~ 530 pc 



( Mizuno et al. 200lh . These circumstances seem 



to point towards a differ ent source of heating. 
Bastian k Goodwin (|2006h reported that, in some 



nearby starburst dwarf galaxies, super star clus- 
ters younger than 60 Myr show a significant devi- 
ation from King and EFF profiles. From the com- 
parison between the observe d luminosity profiles 
and d ynamical simulations, iGoodwin k Bastian 
(2006) proposed that these clusters are far from 
virial equilibrium, because of the likely expulsion 
of the residual gas after the process of SF. We 
therefore suggest, as an alternative to the merger 
hypothesis, that NGC 376 was perturbed by the 
rapid expulsion of gas that followed the initial 
phase of SF. Because of the drop in the initial 
binding energy of NGC 376, a large fraction of 
stars now have velocities higher than the escape 
velocity and are no longer gravitationally bound 
to the cluster. 

We studied how the PDMF varies as a function 
of the distance from the center of NGC 376. We 
found that the PDMF remains flat ((r) = -0.31 ± 
0.27) over the entire c luster, and it b ecomes close 
to the value derived bv lSalpeterl (|l955l ) for the IMF 
in the solar neighborhood only in the clusters's 
tail beyond the tidal radius. If we assume that 
the stars found in the tail of stellar density profile 



13 



Table 2 

Properties of NGC 376, as derived in this paper 



NGC 376 core radius from King's model 
NGC 376 tidal radius from King's model 
NGC 376 core radius from EFF model 
PDMF slope within 20" from the center 
PDMF slope between 40 and 60" 
NGC 376 total mass 
Total mass within 60" 
NGC 376 dynamical time 
NGC 376 escape velocity 
Assumed distance modulus 



r c = 7.95" ± 0.35" ~ 2.34 ± 0.10 pc 
r t = 20.00" ± 0.67" ~ 5.88 ± 0.20 pc 
r c = 8.10" ± 1.58" ~ 2.40 ± 0.46 pc 
r= -0.31 ±0.27 

r = -i.20 ±o.ii 

MTOT rt ~7x 10 4 M Q 
MTOT ~6x 1O 5 M 
T dyn ~ 2.6 x 10 5 yr 
~ 3 km/s 
(m - M) Q = 18.90 



between ~ 20 and ~ 60" belong to NGC 376, then 
almost 90% of total stellar mass of the cluster is 
likely not bound to the cluster. 

Se veral authors (Zhang fc Falll 19991 Chandar et al 



ng 

2006 [ | de Griis fc AndersfeOOftlde Griis fc Goodwin! 
20081: bieles fc Bastianll2008l:IChandar et al.ll2010allbT ) 

have recently analyzed the properties of the young 
clusters (less than few 10 8 Myr) in both the 
Magellanic Clouds with the aim of understand- 
ing the disruption mechanisms of star clusters 
and how long they can survive. However, dif- 
ferent groups derived contradicting conclusions, 
even when using the same data, wit h some fa- 
voring the infant mortal i ty scenario (|Carpenter 
20001 : (Lada fc Lada|[200l IWhitmore et al.l l2007h . 
that predicts that more than 80 — 90% of the 
clusters will be disrupted in a sho rt, nearly mass 



independent interval of time (e.g. IChandar et al 
200fil l2010tJl . and others favo ring a less efficient 



atic study of the dynamical and kinematic prop- 
erties of resolved clusters younger than 10 s Myr 
can provide independent diagnostics to probe the 
time-scale of cluster disruption. 

We thanks the referee for many constructive 
suggestions. We are grateful to J. Anderson for 
the software support, and to N. Panagia, L.R. Be- 
din, and A. Bellini for useful discussions. M.C. 
and M.T. acknowledge financial support through 
contracts ASI-INAF-I/016/07/0 and PPJN-INAF 
2008. 
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